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This paper presents a heavy traffic analysis of the behavior of 
' muhi-class acychc queueing networks in which the customers have 

, deadhnes. We assume the queueing system consists of J stations, and 

• ' there are K different customer classes. Customers from each class ar- 

' rive to the network according to independent renewal processes. The 

, customers from each class are assigned a random deadline drawn from 

a deadline distribution associated with that class and they move from 
station to station according to a fixed acyclic route. The customers 
at a given node are processed according to the earliest-deadline-first 
(EDF) queue discipline. At any time, the customers of each type 
^ ' at each node have a lead time, the time until their deadline lapses. 

We model these lead times as a random counting measure on the 
' real line. Under heavy traffic conditions and suitable scaling, it is 

proved that the measure-valued lead-time process converges to a de- 
terministic function of the workload process. A two-station example 
is worked out in detail, and simulation results are presented to illus- 
, trate the predictive value of the theory. This work is a generalization 

' of Doytchinov, Lehoczky and Shreve [Ann. Appl. Probab. 11 (2001) 

I 332-379], which developed these results for the single queue case. 

S3 ' 

1. Introduction. The last decade has brought dramatic developments in 
communication technology. There are now a wide range of types of com- 
munication services available and an ever increasing demand for those ser- 
vices. An important component of this demand is for real-time applications, 
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that is, applications with specific timing requirements. Examples are video- 
conferencing and video-on-demand, in which the timely delivery of packets 
must be maintained to ensure continuity of the image and sound. Networks 
servicing real-time applications also arise in production systems where the 
orders have due dates or in control systems where there are upper bounds 
on the latency between the occurrence of an event and the control system's 
response to it. 

For queueing networks, the measures of performance are often related 
to system stability, to queue lengths (e.g., the adequacy of buffer space or 
the amount of work- in-process in those buffers) or to waiting times (e.g., 
the delay at a node in the system). For real-time applications, in addition 
to standard measures of stability, queue length and delay, one must also 
be concerned whether individual applications are meeting their timing re- 
quirements. Determining the waiting time distribution is not sufficient to 
determine whether a particular queue scheduling policy can satisfy real- 
time application (customer) timing requirements under various workload 
conditions. 

To model customer timing requirements, we assume that each customer 
arriving to the system has an initial lead time I. If the customer arrives at 
time t, then its deadline is met if it exits the network not later than t + 
To determine whether customers meet their timing requirements, one must 
dynamically keep track of each customer's lead time, where the lead time is 
the time remaining until the deadline elapses, that is, 

lead time = initial lead time — time elapsed since arrival. 

In the study of real-time systems three different types of timing requirements 
are usually distinguished: hard, firm and soft deadlines. Hard deadlines must 
be met or a system failure is considered to occur. These applications arise 
in many control systems, especially in avionics systems or automobile en- 
gine control applications. For a computer system to meet hard deadlines, 
there must be essentially no stochastic aspects associated with the arrival or 
servicing of an application or these quantities must be bounded above. The 
subsequent analysis is based on those bounds with service times taking their 
longest possible value and interarrival times taking their shortest possible 
value. This worst case approach can result in systems functioning at very 
low levels of average case utilization to ensure they can meet application 
deadlines under worst case conditions. 

Applications with firm or soft deadlines are permitted to miss their dead- 
lines, usually with some bound on the size or rate of misses. This class of 
examples include audio and video transmissions, where the end user is able 
to tolerate a small lack of continuity in the sound or image being transmit- 
ted. A firm deadline is one which can be missed and there is no value in 
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completing a task whose deadline has expired, while a soft deadline permits 
lateness and uses the task completed after the deadline. In this paper we 
study the soft deadline case but wish to control the fraction of customers 
that will miss their deadlines and to model the amount of lateness as a 
function of the workload and the scheduling policies used. 

As noted earlier, to study queueing systems in which the customers have 
deadlines, one must attach a lead-time variable to each customer in the sys- 
tem. It is convenient to model the vector of customer lead times at any time t 
as a counting measure on M with a unit atom at the current lead time of each 
customer and total mass equal to the number of customers in the system at 
time t. Exact analysis of such a measure-valued process is intractable; how- 
ever, a heavy traffic analysis can be done. Doytchinov, Lehoczky and Shreve 
(DLS) [5] studied the single queue case in which the customers are processed 
according to the earliest-deadline-first (EDF) or the first-in-first-out (FIFO) 
queue disciplines. DLS proved that under heavy traffic conditions, a suitably 
scaled version of the random lead-time measure converges to a nonrandom 
function of the limit of the scaled workload process, which in the case of 
a GI/G/1 queue, is a reflected Brownian motion with drift. This paper is 
focused on generalizing the results of DLS to acyclic queueing networks. 

There is some other recent work on heavy-traffic approximations for sys- 
tems that handle real-time applications. Van Mieghem [15] studied a single 
server multiclass queueing system. He considered control policies to mini- 
mize the total delay cost incurred by customers over a finite time horizon. 
Markowitz and Wein [12] studied the single machine scheduling problem in 
a manufacturing context that included customer due dates along with other 
model features. Lehoczky [10] gave an informal analysis of a single M/M/1 
queue by constructing the generator for the lead-time vector and showing a 
scaled version converged to a deterministic limit under the EDF or the pro- 
cessor sharing queue discipline. Lehoczky [9] used these results to study the 
performance of a threshold access control policy to reduce customer lateness. 
Lehoczky [11] informally extended the analysis to Jackson networks. 

In this paper we extend the approach and the results of DLS to the case 
of acyclic queueing networks. We assume that a queueing system consists 
of J stations, and there are K different customer classes. Customers from 
each class arrive to the queueing network according to independent renewal 
processes. The customers from each class are assigned a random deadline 
drawn from some deadline distribution associated with that class, then each 
moves from station to station according to a fixed route. The K different 
routes are acyclic, meaning that they visit any of the J stations at most once. 
Upon completion of their route, they exit the system. If the lead time of a 
customer at the time of exit is negative, then that customer is late. While 
lateness is permitted, we wish to determine (and ultimately to control) the 
fraction of customers that will exit late. Each station will process some 
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subset of the K different customer classes. At each station, customers are 
queued in lead-time order and preemption (preempt-resume) is allowed. We 
assume there is no overhead associated with preemption. In the special case 
that all customers are assigned zero initial lead time, the EDF policy used 
in this paper becomes first-in-system- first-out (FISFO). Our analysis thus 
provides information about the time customers spend at various stations 
and in the system when FISFO is used. 

We study the behavior of these acyclic networks under heavy traffic 
conditions. That is, we consider a sequence of queueing systems in which 
the traffic intensities at each node approach 1. We prove that if the suit- 
ably scaled J-dimensional workload process converges appropriately (see 
Assumption 2.1 and following discussion), then a suitably scaled version of 
the ii'-dimensional lead-time measure process converges to a deterministic 
function of the workload. 

This paper is organized as follows. Section 2 presents the model, notation 
and assumptions. Section 3 presents the measure- valued processes associated 
with customer lead times and the concept of the frontier processes. This 
section also states the relationship between the limiting workloads and the 
limiting frontiers, a major result of this paper, which is proved in Section 
6. The frontier process was defined in DLS for the single-queue case and 
formed the basic methodology used to analyze EDF and FIFO queues. This 
methodology is generalized to the network case in Section 4, which provides 
several technical results needed for the analysis. Section 5 shows how the 
equations which provide limiting workloads as a function of limiting frontiers 
can be inverted, so that one can determine the frontiers as a function of the 
workloads. Section 7 presents a simple but interesting worked-out example. 
It also presents simulations which illustrate the accuracy of the heavy traffic 
theory. 

2. The model. 

2.1. System topology. We consider a sequence of queueing systems in- 
dexed by n. It is assumed that each system is composed of J stations, 
indexed by 1 through J, and K customer classes, indexed by 1 through 
K. Each customer class has a fixed route through the network of stations. 
Customers in class k, k = 1, . . . , K , arrive to the system according to a re- 
newal process, independently of the arrivals of the other customer classes. 
These customers move through the network, never visiting a station more 
than once, until they eventually exit the system. However, different customer 
classes may visit stations in different orders; the system is not necessarily 
"feed-forward." We define the path of class k customers as the sequence of 
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servers they encounter along their way through the network and denote it 
by 

(2.1) Vik) = {jk,i,jk,2, ■ ■ ■ JkMk))- 

In particular, class k customers enter the system at station jk^i and leave 
it through station jk,m{k)- If j is a member of the list of station indices in 
V{k), we shall write j € 'P{k). 
For j = 1, . . . , J, we define 

(2.2) C{j) = {Indices of customer classes that visit station j}, 

(2 3) IC ( ') — I ^^^^^^^ customer classes that enter 1 
^ ' ^ [ station J from outside the system J' 

We assume that for every j S {1, . . . , J}, C{j) ^ 0. For j = 1, . . . , J and 
k GC{j), we further define 

{Indices of stations visited by customer class 
k before visiting station j, or if customer > . 
class k enters the system at station j ) 

We assume that the network is connected, in the sense that for any two 
stations, there is a way of reaching one station from the other by following 
fragments of paths of the form (2.1), not necessarily in the forward direction. 
The network topology captured by (2.1)-(2.4) does not depend on n. 

2.2. Arrival times, service times and lead times. The customer inter- 
arrival times are a sequence of strictly positive, independent and identically 
distributed random variables Ujl^"'\ i = 1,2,..., where the subscript k in- 
dicates the customer class and the superscript (n) indexes the queueing 
system. These are independent across k as well as i. We assume that 

(2.5) Al-)4(Enl'("V\ 4")^(Varul;("V/^ 

are both defined and finite. 

The customer service times are a sequence of strictly positive, indepen- 
dent and identically distributed random variables i"]' ^ = 1,2,..., where j 
indicates the station of service and k £C{j) indicates the class of customer. 
The superscript i indexes the order of arrival of customers of class k to the 
system, which may be different from the order of arrival of class k customers 
to station j. The random variables fc"] are independent across k and j as 



well as i, and they are independent of the inter- arrival times. We denote by 
Vu^"'^ the service times of customers of class k at station j, with index i indi- 

eating the order of arrival to station j. The random variables {v^^j ^}^]^ are 
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thus a random permutation of the random variables {vq ^j}^]^. Under the 
EDF service disciphne described below, the index i of arrival at station j of 
a customer of class k is independent of f^' j^''. Therefore, v'^^^\ i = 1,2, ... , 
are also independent and identically distributed, with the same distribution 
as the random variables t'oi"]- ■'^^'^ j = ■ ■ i J ^^'^ k £C{j), we assume that 

^ /?e = (Var.g)^/^ = (Var<("))V^ 

are both defined and finite. 

Each customer in class k arrives with an initial lead time L^^^"^^ having 
distribution 

(2.7) P{4'(")< V^y} = Gfe(y), 

where is a cumulative distribution function. These lead-times are inde- 
pendent and identically distributed across i. They are independent across 
k and independent of the interarrival and service times. Note that Gk does 
not depend on n; the lead-time distributions dilate at rate y/n as oo. 
We assume that for k = 1, . . . ,K , 

(2.8) y^^sup{yGM;Gfc(y)<l}<cx). 

We also assume that for every n, the sequences {u^'^"^}, {foi"]} ^'■^'-^ {-^l'^"''} 
are mutually independent over j E {1, . . . , J}, k G C(j) and i = 1, 2, 

2.3. EDF discipline. Customers are served at each station according to 
the EDF discipline. That is, the customer with the shortest remaining lead 
time, regardless of class, is selected for service. We give the analysis for the 
case of no preemption. If preemption is permitted and we assume preempt- 
resume, then an obvious simplification of the analysis gives the same results. 
There is no set-up, switch-over or other type of overhead. Late customers 
(those with negative lead times) stay in the system until served to comple- 
tion. We assume the system is empty at time zero. 

2.4. Unsealed queueing processes. For each station j = 1, . . . , J and cus- 
tomer class k = 1, . . . ,K , we define 

m 

(2.9)5r'('^)^E%'^"^ 

i=l 

= Time of arrival to the system of the mth customer of class k, 

(2.10|(")(t)^max{m;5r'^'^)<t} 

= Number of class k arrivals to the system by time t, 
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(2.11^["j(t) = Number of class k arrivals to station j by time t, 
(2.12^^"j(t) = Number of class k customers at station j at time t, 

= Number of customers at station j at time t, 

i=l 

= Work for station j associated with the first [t\ customers 
of class k to arrive to the system, 



i=l 



= Work associated with the first \t\ customers of class k 
to arrive at station j. 

(Here and elsewhere we use the convention X]i=i = J2k&0 = 0-) We further 
define the netput at station j to be 

(2.16) ivf)(t)^ E vi-^iAt']{t))-t, 

kdCii) 

and the cumulative idleness at station j to be Ij^\t) = — info<s<f ivj"'^(s), 
which is nonnegative because Nj{0) = 0. Finally, the workload at station 
j is the amount of time it would take to serve all customers at station j 
to completion, assuming no new customers arrive, and this is W^''\t) = 

ivf)(t) + /f^(t). 

2.5. Scaled queueing processes. In order to obtain a limit as n ^ c>o, it 
is necessary to scale and sometimes center the above processes. We define 

= ^iA^k\nt) - Ai")nt], 4 -^[A^M " X^k^nt], 
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2.6. Heavy traffic and convergence assumptions. For j = 1, . . . , J and k G 
C{j), the traffic intensity of class k customers at station j is p^^j = A^^V/^i";-) 

and the traffic intensity at station j is p^"'^ = J2keC(j) P^kj- assumed that 
for all j, 

(2.17) 7,-^ lim ^^(1-^;.")) 

exists. Furthermore, it is assumed that for ah k and j satisfying k G C(j), 



(2.18) 



Afc = lim A^"\ pk,j = 

(n) 

n— +00 n— >oo 



= lim 



lim /il"], 
hm /3["^ 

I— >oo '^'J 



are all defined and and pf^j are positive. We define the limiting traffic 
intensities pkj = Xk/pkj and pj = J2kec{j) Pk,j- We impose the usual Lin- 
deberg condition on the inter-arrival and service times: for j = 1, . . . , J and 

k£Cij), 



lim E 

n— >oo 



(2.19) 



: lim E 

n—>oo 



-1\2 



{|«^(")-(A^')-i|>Cv/ii} 



-1\2 



^{i<^'-(A'i::;)-M>c^} 







VoO. 



In what follows, the symbol =^ denotes weak convergence of measures on 
the space Ds[0, oo) of right-continuous functions with left limits from [0, oo) 
to a Polish space 5. The topology on this space is a generalization of the 
topology introduced by Skorokhod for D5[0, 1]. See [2] for details. We take 
S = M (or W^, with appropriate dimension d, for vector-valued functions) 
unless explicitly stated otherwise. 

Theorem 3.1 [14], together with (2.19) and the independence assumptions 
of Section 2.2, implies that for every j = 1,. . . ,J and k G C{j) and every 
y <V*k-, we have 



(2.20) rg/t;?/) 



1 



[nt\ r 

E 

1=1 L 



„,*.(")l[ 

^0,fc,i^{L^.(")<^y} 



1 



(n) 



Gk{y) 



TlAt;y), 



where Tj^ j{t;y) is continuous in t. Putting y = yl into (2.20) and using the 



fact that the sequences {foi"]}i^i {^A: ^'''ii^i have the same distribution, 
we conclude that for j = 1, . . . , J and k £ C{j), 



(2.21) 



Vr 



(n) 



0,k,j 



k,ji 



(n) 



where ■ is a continuous process. In fact, if /J^j > 0, then ^ is a Brownian 
motion. Similarly, Theorem 3.1 of [14] and Theorem 14.6 of [2] imply that. 
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for every k, there exists a continuous process such that 

(2.22) ^fU^^ 

By (2.21), (2.22) and a standard argument (see, e.g., [8], Corohary 3.2), we 
have 

(2.23) M^i^it) ^ j^[vi:^,iAt'\nt)) - npgt] ^ iV4*^.(t), 



where j is continuous. We also make the following convergence assump- 
tion: 

Assumption 2.1. For every j and k £ C{j), there exists a continuous 
process A^. ^ such that 

(2-24) ^Ij^^lj- 

There exists a J-dimensional continuous process (W^ , WJ, . . . , Wj) such that 
(2.25) iwl''\wi''\. . . , W}""^) iWl,W^, W*j). 

In feed-forward networks, (2.24) and (2.25) hold under FIFO (see [13]) 
and EDF (see [17]). Because our network is not of the feed-forward type, 
there are no known general conditions which guarantee (2.24) and (2.25). 
However, the literature contains a number of special cases of our model in 
which (2.24) and (2.25) hold. Rather than take one of these special cases as 
a starting point, we choose to begin with Assumption 2.1 because this is all 
we shall need in order to obtain convergence of scaled lead-time profiles. 

3. Measure-valued processes and frontiers. Our goal is to obtain a char- 
acterization of the lead-time profiles of the customers queued at the J sta- 
tions in the system in terms of the limiting workload process {W^ , W2 , • • • , Wj) 
in (2.25). These lead-time profile processes are measure- valued. More pre- 
cisely, they take values in the space M of finite, nonnegative measures on 
B{M), the Borel <T-algebra on M, equipped with the weak topology. In what 
follows, we shall denote by A^"^ the J- fold product of Ai (with the product 
topology). For a Borel set B cM, we set 

{Work for station i represented by 
class k customers at station j with > , 
lead times in B at time t J 

f3 2) W'^'^^^ (t)(B) — /^°^^ station i represented by customers 1 
^ V 7V 7 |g^^ station j with lead times in B at time t J ' 
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Then wj^j {t){B) = yvl'^j^\t){B) is the work at station j represented by class 

k customers at that station with lead times in B at time t, and Wj^\t){B) = 

is the work at station j corresponding to all customers at that 
station with lead times in B at time t. We also define 

(3 3) Q*'"^(i)(-B) ~ [Number of customers at station j 

\ with lead times in B at time t 



{Number of class k customers arriving ' 
to the system by time t and having 
lead times at time t in B, whether or 
not still in the system at time t 



and 



{Work for station j associated with customers 
of type k arriving to the system by time t and 
having lead times at time t in B, whether or 
not still in the system at time t 

The scaled versions of these processes are 

Wif\tm^^Wif\nt){V^B), 



n 



y^l1{tm^^w'i']{nt){v^B) 



Qf(t){B)^^Qf{nt){^B\ 



At\tm^^At\nt)[V^B\ 



in 

We introduce frontier processes 

{Largest lead time of any class k customer 
who has ever been in service at station j, 
or \^y1 — t if no such customer exists or 
if this quantity is larger than the former one , 

^ ' ' kadi) "'^ 
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The scaled versions of these processes are 

(3.8) = ^i^SM, Fj^'Ht) = ;^^^"^(-*)- 

The next step is to define a set D which contains the J-dimensional 
vector-valued process To do this, we begin with a 

permutation vr = (vri, 7r2, . . . , vrj) of the integers (1, 2, . . . , J). Given such a 
permutation and an integer m G {1, . . . , J}, we define 

(3.9) ICl,,^,{j) ^{k€ CU);S{k\j) C {vri, . . . ,7r^_i}}, 

(3.10) J--„i ^ {j; /C- _i(i) / 0} \ {vri, . . . , vr^.i}. 

By convention, if m = 1, then {vri, • • • , T^m-i} = 0, ^oU) — ^o(i)) the set of 
indices of customer classes which enter the system at station j , and = , 
the set of stations which serve as the entry point for at least one external 
arrival process. These two sets do not depend on the permutation vr. Subse- 
quent sets do. The set lC'^_i{j) is the set of all customer classes that visit 
station j and visit only stations in the set {vri, . . . , vTm-i} before arriving 
at station j. We say it is the set of customer classes which reach station j 
through {vri, . . . , VTm-i}. The set iZ^-i is the set of all stations j not in the set 
{tti, . . . ,TTm-i} which are visited by at least one customer class of the type 
just described. We say that is the set of stations which can he reached 

through {vri, . . . , TTm-i}. Note that both /Cm_i(j) and J^-i depend only on 
(vTi, . . . , TTm-i), not the full permutation. Thus, we shall sometimes write 

1C^Z-r^"'^^\j) and J-ifii-'^'"^'^ instead oi }Cl^^^{j) and J^^^. Finally, we 
define 

, ^^-^ n = {vr; vr is a permutation of 1, . . . , J 

and TTm G Jm-i all m = 1, . . . , J}. 

In other words, 11 is the set of all permutations vr = (vri, • • ■ > ^rj) such that, for 
each m, the station vr^ can be reached through {vri, • • • >7rm-i}- For vr G 11, 
we set 

(3.12) D'' = <^ y G ]R-^;y^j >•••> and < max Vm 

(3.13) D^[j D"". 

Lemma 3.1. For all t> id, the random vector {F[''\t),F^''\t), . . .,Ff\t)) 
takes values in the set D. 

Proof. We must construct a permutation vr = (vri, . . . , vrj) G n such 
that Fif{t) > Fif{t) >■■■> M'iHt) and F^^(i) < maxfceyc^_,(.„) for 
every m. We do this by induction. 
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We note first that because customer class k visits consecutive stations 

(3.14) V^yl > Fj^lint) > F£.) ^^(nO > • • • > (nt). 

This imphes that the largest frontier must be at a station which is in Jq, 
the set of stations that have arrivals from outside the system. We select a 
station tti G jTo whose frontier Fn"^^ [nt) is maximal. If this maximal frontier 
is the lead time of a customer which has been in service, we may choose vri 
to be the station where that customer entered the system. If the maximal 
frontier is not the lead time of a customer who has been in service, then it 
is of the form \/ny^ ~ ''^^ fo'^ some customer class k. In this case, we choose 
TTi to be the station where this customer class enters the system, so that k 
is in A^o(7i"i)> the set of customer classes that enter the system at station tti. 
In either case, we obtain 

(3.15) max ^/nyl> F^l\nt) >xaayiF'f'\nt). 
fceKo(i"i) jyTi 

For the induction hypothesis, we assume for some m G {2, . . . , J} that we 
have constructed vri, . . . , vTm-i such that: 

(i) for each i <m — 1, station tTj is reached through {vri, • ■ • , T^i-i}, 

(ii) for each i <m— 1 , 

(3.16) max V^yl > F^f (nt) , 

(iii) we have 

(3.17) FW(nt) > • • • > Fi±^{nt) > max &\nt). 

J^|7ri,...,7r™,_i} 

If the maximal frontier among Ff'\nt) for j ^ {vri,... ,7rm_i} is the lead 
time of a customer that has been in service, we may choose vr^ to be the 
station where that customer first reaches a station not in {vri, • • • iT^m-i}- If 
this maximal frontier is not the lead time of a customer who has been in 
service, then it is of the form y/ny^ — nt for some customer class k. In this 
case, we choose vr^ to be the first station not in {vri, • • • iT^m-i} reached by 
this customer class. In either case, we obtain 

(3.18) max V^yl> Fi'l^{nt)> max Ff{nt). 

The induction step is complete. 

Once the induction has concluded, we have constructed a permutation 
satisfying properties (i)-(iii) for m = J + 1. Dividing (3.16) and (3.17) by 
y/n, we obtain the desired properties for the scaled frontiers. □ 
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For k = 1, . . . , K , we define 

(3.19) Hk{y)= {l-Gk{x))dx, y G M. 

This function is strictly decreasing on (— cx3,y^J, mapping tliis half- line onto 
[0,cx)). We next define $ = ($i, . . . , $j) : M-^ ^ [0,oo)-^ by 



(3.20) 



^j(yi,---,yj) - Pk,j 



Hk{yj)-Hk\ mill 

»6<S(fc|j) 



J = 1,...,J. 



In the above definition and in all that follows, the minimum taken over the 
empty set should be interpreted as oo. The main results of this paper are 
the following two theorems. 



Theorem 3.2 (Convergence of scaled frontiers). The function ^ is a 
homeomorphisni of D onto [0,00)''. With 



(3.21) 
we have 
(3.22) 



{F^,...,F})^'^~\W^,...,W: 



(Ff 



(n) p(n)x 



We define W*it) = (Wi*(t), . . . , W}(t)) and Q*(t) = (QI (t), . . . , Q}(t)), 
which take values in A^'^, the set of J-dimensional vectors of measures on 
M, by specifying their values on half-lines of the form (y,oo) for all y E M. 
This is done for j = 1, . . . , J by the formulas 



(3.23) W*{t){y,oc)^ ^ p^j 



Hk{yVF;{t))-Hk( mm F*{t) 
teS(k\j) 



(3.24) Q*(t)(y,oo)^ 

fc6C(j) 



Hk{yVF*{t))-HJ mm F*{t) 
■' \*e'S(fc|j) 



Theorem 3.3 (Convergence of scaled workloads and queue lengths). 
We have 

(3.25) (w/"\...,>vj"))^>V*, (qS"),...,q5^))^Q*. 



The weak convergence in (3.25) takes place in Dj^j[Q, 



00]. 
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4. Customers behind the frontiers. In this section we prove the crucial 
observation that both the number of customers at each station j with lead 
times smaller than or equal to the current frontier Fj"'^ (t) and the work for 
the system associated with these customers are negligible. This is done in 
several steps, leading to Corollary 4.7. Along the way, we show tightness of 
the scaled frontier processes (Lemma 4.6). Both these results will be used 
in Section 6 to prove Theorem 3.2. 

Proposition 4.1. Let j = 1, . . . , J, k e C{j), -oo <yo<yl and T >0 
be given. As n ^ oo , 

(4.1) sup sup \V^''l.{t){y,oo) + pkj[Hk{y + V^t)-Hk{y)]\^0, 

yo<y<ylo<t<T ' " 

(4.2) sup sup \Al!'\t){y,^) + Xk[Hk{y + VEt)-Hk{y)]\-^0. 

yo<y<yl o<t<T 

The processes in Proposition 4.1 do not take departures into account. 
Because this proposition is concerned only with arrivals, its proof can be 
given following the proof of Proposition 3.4 of [5]. We do not repeat that 
proof here, but instead give a heuristic argument. Let us first consider (4.2), 
which asserts that asymptotically, the "density" of the measure- valued pro- 
cess ^f)(t) is the same as the density 

(4.3) xthH'.iy + V^t) - H'M] = >^t\Gk{y + V^t) - G,{y)]- 

In order for a class k customer to have lead time y at time t, it must arrive 
at some time t — s prior to t and be assigned lead time y + s upon arrival. 
If Gk has a density, then the density of the assigned lead-time distribution 
is ;^G'^(^^) [see (2.7)], and multiplying by the arrival rate A^"'^ we obtain 
the density of class k customers with lead times y: 



n Jo \ y/n 



n J \ \ n 



The heavy traffic scaling considers the density of times the actual number 
of customers whose lead times are y = at scaled time t = This density 

is the right-hand side of (4.3). (The Jacobian ^ = ^/n is canceled when 
we divide the customer count by \/n.) Under the heavy traffic scaling, the 
work brought by customers of class k to station j is the average work per 
customer, (^fc"])"^, times the number of customers. Multiplying the right- 
hand side of (4.3) by we obtain p^"] [^^(t/ + ^t) - Gfc(y)], which 
explains (4.1). 
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Corollary 4.2. Let j = 1, . . . ,J , ke C{j), -oo < yo < Vk "-"^d T>0 be 
given. As n oo, 

sup sup v^"2j(0{y} -^o> 

yo<y<yto<t<T 

(4.4) yu-y-yfe _ _ ^ 

sup sup 4"^(t){y}^0. 

yo<y<ylo<t<T 

This coroUary is a consequence of the fact that the hmiting measures for 
^o^kji^) •^k^\^) have densities and, hence, do not charge points. Its 
proof is similar to the proof of Corollary 3.5 of [5], and we refer the reader 
there for details. 

Lemma 4.3 and Corollary 4.4 generalize Proposition 3.6 in [5] to the case 
of acyclic networks. We provide the details of these proofs. 

Lemma 4.3. For all k and j G 'P{k), we have 

(4.5) H-oo,Fi5) ^0, W,(")(-oo,Fi:;.)) ^0, 
and for every station £ satisfying j G S{k\£), 

(4.6) n?'("\-oo,FiJ)^0. 

The weak convergences in (4.5) and (4.6) take place in L'k[0,oo). 

Proof. We fix ko and prove (4.5) and (4.6) by induction along V{kQ). 
Let jo G V{ko) be given and make the induction hypothesis that (4.5) and 

(4.6) hold for k = ko and every j £ 5(A;o|jo). If jo is the first station in Viko), 
this hypothesis is vacuous. 

From the induction hypothesis (4.6) we have 

for every j G S{ko\jo). For such a station j, we have -F)^" j- > -^fc" io hence, 

(4.7) H5(")(-oo,i^i)^0. 

We now prove (4.5) and (4.6) for k = kQ and j = jo- Toward this end, we 
define 

4:i{t) = sup{. < t- Q'^\ns){-^,Fl^l{ns)) = 0} 

= sup{. < t; q(.;;)(.)(-oo,f£)^(.)) = 0}, 

and note that 
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Indeed, from the definition of ''"i" 

(4-9) 25?(<kw-)(--'<L(-So(*)-))=o 

and the only way in which F^-^ can jump up at time 't-t^^ is that a 
customer of class with lead time greater than any customer of this class 
who has ever been in service at station jo begins to receive service at time 
nT^-^{t). But then, by the EDF service discipline, 

2i?(<kw)(--'<L(-i:i(^)))=o 

and, consequently, the value of Qjo^(")(~°0' -^i^jo^")) remains zero at least 
until the next arrival to server jo after time ^Tk^^jg{t), which contradicts the 
definition of T^j^^it)- Thus, F^]j^^ cannot jump up at time nrj^^j^it), so (4.8) 
follows from (4.9). For T^"jg(0 <s<t, we have 

(4-10) Fi^lins) = Ft:l^inri:l{t)) - n(. - rj^^^it)). 

The first step is to prove 

(4.11) *-^£oW^O' 
and, subsequently, to upgrade this convergence to 

(4.12) V^(t-ri;;)^(t))^0. 

The convergences in (4.11) and (4.12) are in Z)ig[0,oo). The key inequality 
is 

w];)(nt)(-oo,F£)^(ni)) 

< if (nt) + J2 '^t^ + ^^"^ (^*) 

(4.13) i65(fco|jo) 

+ i^f^(nt)-n(t-r£)^(t))+i?W(nt), 



where the terms on the right-hand side are defined below. The first term, 

H^-\nt) 4 wj;)(nr(;:).^(t))(-oo,F£).^(nrj^^^ 



accounts for the work arriving to station jo by time '^''"i" jg (^) • A typical 
summand in the second term. 
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is the work for station jo at upstream station j which is ahead of F^^^j^ (""^fco jo ^^-^ 

at time f^'r'^j^ and, hence, has the potential to arrive at station jo by time 

nt ahead of F^^-^{nt). If jo is the first station in V{ko), these terms do not 
appear. The third term, 

oo 

D^'^\nt) =51 VfcojO^{nr<"'. (t)<S?'(")<ni}^{L!.'("'-(ni-5!''"')<F("', (nf)} ' 

is the work of type ko arriving to the system during the time interval 
(nr^"^j^ , ni] with lead time upon arrival that puts it ahead of Fkgjg{nt) at 
time nt. A typical summand in the fourth term. 



oo 



(nT^^j^,nt]. The fifth term, —n(t — rj^^j^it)) , is the work accomplished by 
the server during the time interval (nr^"^^-^(t), nt], some of which may be 
devoted to a customer already in service at time t^t^^I- (t) whose lead time 



is the work of type k ^ arriving to station jo during the time interval 
^,jo'^*]- '^^^ ^^^^ ~^(*~^£jo( 

is greater than or equal to F^^-^^{nT^^-^^{t)) and the remainder of which at 
each time s G {'^koja^A ™ust be devoted to customers with lead times less 
than f["'\ (ns) . The final term, 

i?(")(nt) ^ maxl^J); 1 < i < ^'J^lnt), A; G C(jo)}, 

is an upper bound on the amount of work that can be devoted to a customer 
already in service at time '^''"i" jqI^)- If there is preemption, this final term 
does not appear. 

We fix r > and estimate the terms appearing on the right-hand side of 
(4.13). The terms o(^/n), 0{y/n) and 0(n^^/^) in the following argument 
depend on T but not on t G [0,r]. For t £ [0,T], (4.8) implies 



(4 14) ^-""-^ 



because in (2.25) is continuous. Furthermore, 



)(nt) < mc«c^n^;(")(n.)(-oo,F£)^(ns)) 
^^■^^^ Wax^)^f )(.)(-oo,f£^(.)) = o{V^) 
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because of (4.7). The last term, R^'^\nt), satisfies 

(4.16) i?(")(nt) < max [wl^'^ns) - wl''\ns-)] = oiV^) 



as in (4.14). We also have 



Jn) _ 1 \ _ ^ IJn) 



(n) ) I ""fejo (n) 

^S(^41(nt))-i>S(^41(<lw) 

i>S(;^^H(*) + ^*') 



+ i^If [4S(*) - ^So(A:im + Piln(t - r£)„(t)). 
(4.17) 

From (2.21) and (2.24), we obtain 

(4.18) Kl-\nt) = pg^n(t - r£)^(t)) + O(V^). 

Finally, we estimate D^''\nt). For this we choose y < yl and divide the 
analysis into the two cases 

(4-19) t-r£^^(t)<^(yl.„-y) 



<nt. 



and the complementary case 

(4.20) „^W.^(f) + ^(y*^_y) 

In the former case, t — Tko\o(^) ~ 0(n^^/^). We show this is also true in the 
latter case. Under condition (4.20), because 

4:iint) = <)„(<)„ (t)) - nit - ri:lit)) < V^y^ - n{t - r£„(0), 
we have 

D(")(ni) 



oo 

<: \^ IT 

- Z^^O,fcojo^{nr<"', (t)<5'''"^<nr'"'. (i)+v^(y* -y)} 
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(4.21) 



(n) 



In) 



+ E 



i=l 



E 



i,(n) 



(n) 



GkM 



+ - 4?("-So w + ^(^^0 - 2/))] 



-in) 



1 



p(n) 



0,fco JO 



+ 



(n) 



(n) Jn) 



n 



From (2.20), (2.22) and (2.23), we obtain 

(4 22) ^^"^(-*)^"G.„(,)^i^^.„(t-r£Ui)) 

+ Mi-GM)p'-ZMo-y) + o{V^)- 

Substitution of (4.14)-(4.16), (4.18) and (4.22) into (4.13) yields 



0<w]:)(nt)(-oo,F£)./nt)) 



(4.23) <n(t-r£).^(t)) 



E pBo-1-(1-G'.o(2/))p£!,o 



(n) 



< 



LfcecOo) 



+ o(^/^) 



(n) 



where the last inequality follows from (2.17). Assume for the moment that 

y < yl^. Then (1 — Gka{y))p^^^ ■ is strictly positive and bounded away from 
zero uniformly in n. This implies 
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Hence, (4.11) holds. 

Armed with (4.11), we return to the weaker assumption y < y^^ and use 
the differencing theorem in (4.17) (see, e.g., Theorem A. 3 of [5]), (2.21) and 
(2.24) to obtain 

(4.18') Kt\nt) = p^:}n{t - rQ^^{t)) + o{V^) 

in place of (4.18). Similarly, (4.22) becomes 



(4.22') 



Z)W(nO<nG,„(y)pg^.„(t-ri").^(t)) 



If (4.20) holds, we may substitute (4.14)~(4.16), (4.18') and (4.22') into (4.13) 
to obtain 

0<w];)(nt)(-oo,Ffc„,,„(nt)) 

(4.24) < -n{t - r£) „ (t)) (l - (,)p£,„ - pgo) 

+ V^{i - GkM)ptL{yl -y) + o{V^). 

Assume for the moment that y < y^^. Then (4.24) implies 

(4.25) V^{t - r£).„(t)) < Cniyl -y) + 
where 



JO 



The constants C„ converge to a finite limit as n ^ oo, which is bounded 
uniformly in y < y^^. Since y may be arbitrarily close to y^^, by (4.25), if 
(4.20) holds, then (4.12) holds also. Similarly, if (4.19) holds, then we have 
(4.25) with C„ = 1 and again (4.12) follows. 

From (4.12) we immediately obtain (4.5) for k = ko and j = jo, because 

Q^;:\nt){-oo,Fi:l{nt)) 

We divide this by y/n and use the differencing theorem, (2.24) and (4.12) to 
obtain the first relation in (4.5). For the second part of (4.5), we set y = y^^, 
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so that case (4.19) is vacuous. Then (4.20) implies (4.24), which we divide 
by ^Jn. The conclusion follows. 

It remains to prove (4.6) for k = /cq, j = jo and i satisfying jo G S{kQ\tj. 
From (4.8) and (4.10) we see that all work at station jo for station I present 
and having lead time in (— oo, -F^"j^(nt)) at time nt must arrive in the time 
interval [nT^^-^^{t),nt\. It follows that 

(4.26) Wf'(")(nt)(-oo,F£)^(nt))<i7f)(nt)+ 

feeC(jo)ncW 



where 



[nt) ^ (nrj:^^^ (t) ) (-oo, Fj£^^ (nri^l (*) ) ) 



accounts for station £ work arriving to station jo at time 'nT^j^{t), and 

oo 

T^i"^) A „/.io, Wit 



is the class A;, station I work arriving to station jo during the time inter- 
val {nT^-^{t),nt\. The random variables v^^''^'^'^\ i = 1,2, ... , are a random 



permutation of the class k, station I service time random variables Vq'^^\, 
i = 1,2, . . . . The latter are indexed in order of arrival to the system; the 
former are indexed in order of arrival to station jo. Because the index i of 
arrival of a customer of class k to station jo is independent of the service time 
"^fc^"'^"^ of that customer, the sequence v]^'^£'^"'\ i = 1,2, . . . , is independent 
and identically distributed, with the same distribution as v^^^l, i = 1,2, . . . . 

We bound the terms appearing on the right-hand side of (4.26). We have 
first of all that 

1 .An) 



HPint) 
< > max —=Vr.\f 

ss^^„5^.(^;:,(i4"'(™))-t>„<£(i4"'(».-)).^; 

But the process 
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converges weakly in L>]8[0,cxo) to the continuous process V^ii^ks) [see (2.21) 
and (2.22)] and, thus, its maximum jump over s G [0,T] converges to zero. 
It follows that 

(4.27) max -^Hj;''\nt) => 0. 

^ ' o<t<T y/n 



Let us now define 

[nij 



^ . , , 1 L"*J / 



k,l („) I' 



which satisfies V"/"/"^ ^ 14*^, where V^*^ is a continuous process [cf. (2.21)]. 



Then 



The right-hand side converges weakly to zero in -Dir[0, oo) because of (4.11), 
(4.12), the continuity of ^ and Assumption 2.1. In particular, 

(4.28) max -^K^^jint) ^0 VA: G C(jo) nC(^). 

o<t<T -/n i^^*-^ ' 

From (4.26)-(4.28), we have 
max^>^"^'^^W(-oo,Fi:i(t)) = ma^^^ 

This gives us (4.6) for k = k^, j = Jq and £ satisfying jo € 5(fco|^). □ 

Corollary 4.4 (Crushing). For every j , we have 

(4.29) Q(.")(-cx),Ff^(t))^0, >v/"Voo,Ff^(i))^0. 
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Proof. From (4.5) we have 

Q^;\-oo,Fj-\t)) = max Q^\- oc, Fj^;^ (t)) 0, 

^(-cx),Ff = max w\'^\-oo,fI^] {t)) ^ 

■' ■' k£C{j) " 



□ 



Corollary 4.5. For all stations j and i for which there exists a cus- 
tomer class k such that j £ S{k\£), we have 

Proof. It suffices to show that for every k E C{j) nC{i), we have 

For such a k, let J^"-(nt) denote the indices i, according to the order of 
arrival to the system, of the class k customers at station j at time nt with 

j 



lead times in {—oo,F^^\nt)). Then 



n ^ o,k,i- 



Let \J)!''j{nt)\ denote the cardinality of For each positive inte- 

ger m, let = P{|j[."j(nt)| = m}. Let Fq^^ be the zero measure, and for 

(n) 

each positive integer m, let F'rn' denote the measure induced on M by the ran- 
dom variable ^ Y]. , fni"]. Conditioned on |jl'^](nt)| = m, the distribu- 



tion of ^^'^ same as the distribution of Ym!=\ ^ofc'l- Because 

the random variables fQ^;?] do not enter the determination of the indices 

which below to J^, the measure induced on M by {\i^i^\nt)\>\} ^ieaj^"'{ni) ' 

is Y^^=id'^mF'^m ■ The set of probability measure {P-m^} is tight because for 
m > 1 and K > (/i^"])^^, we have 



1=1 >i 1=1 



1 m 



< ri IE 



m . 
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and this can be made arbitrarily small, uniformly in m and n, by the choice 

of K. Consequently, the set of probability measures {Z)m=o '^"*^^m''}n>i 
also tight. In particular, given e > 0, there exists K > such that 

for all n > 1. 

According to the first part of (4.29), (nt)| =^ 0, and hence there is 

an integer N such that 



for all n> N. Therefore, for n> N, 



nwlf\t){-oc,Fj^\t))<e} 

1 - ■s-^ j,(n) 1 |r,(n). 



> 



> 1 - 2e. 
This establishes the corollary. □ 

The following lemma gives a tightness bound for the scaled frontiers. 

Lemma 4.6. For every T >0, e > 0, j G {1, . . . , J} and k G C{j), there 
exists y G (— oo,y^) such that for all n, 

(4-30) p{^mf^Fi;;)(t)<y}<e. 

Proof. As in the proof of Lemma 4.3, we fix k and proceed by induction 
along P(fc). Let £ G T^ik) be given and assume that for j G S{k\i),T > 
and e > 0, the corresponding y satisfying (4.30) for all n can be found. In 
particular, no assumption is necessary to analyze the first station in V{k). 

We first argue that 

(4.31) wi^-\t) > vg,,(t) (^irA*)'^.,^(\^,)^S + '^(1)- 

Indeed, the workload at station i is at least as great as the workload brought 
to station £ by class k customers with lead times in (F^"^(nt), minjg5(^.|£) {nt)). 
None of class k customers who have arrived to the system by time nt and 
have lead times at this time greater than F^^{nt) has ever been in service 
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at station i by time nt. Thus, every such customer is either in queue at 
station i or in queue at some station jo G S{k\i). By Lemma 4.3, for such Jq 
we have 

)^°'(")(t)(-oo,^.mm Fg(i))<^^^^ 

so the difference between Vq^^ g{t){F^^^ {t) ,minj^g(^^g-^ F^"'j{t)) and the scaled 
workload for station £ associated with class k customers already present at 
£ with lead times in {pj^"'^ {nt) , mm pj^^'j {nt)) is of the order o(l). This 
justifies (4.31). 

Fix r > and e > 0. By the induction hypothesis, there exists yi < 
such that, for all n, ¥'{An) > 1 — |, where 

An = { inf min F^^ht)>yi\. 

By (2.25) and the continuous mapping theorem, 

sup W^/")(t)^ sup W^{t). 

0<t<T 0<t<T 

Inequality (4.31) and the fact that limj^-^_oo Hi^{y) = oo enables us to choose 
2/2 < 2/1 such that, for all n, F{Bn) > 1 — |, where 

Bn^{ sup vtVt)fe:;(t),.min Fg^(t))< 
On An n Bn, for < t < T, we have 
(4.32) sup vS,(t)(Fg(t),yi] < VHM. 

0<t<T 

By Proposition 4.1, we can find N such that for all n> N, F{Cn) > 1 — |, 
where 



Cn = I sup sup Vo"fc + Pk/Wkiy + Vni) - Hk{y)] 

^y2<y<yl0<t<T ' ' 



For n>N, P(A„ n B„ n C„) > 1 - f . By (4.32), on A„ n S„ n Cn we have 



^sup^vJJf(i)(y2,yi]I|j?(n)(t)<y2} < ^Hk{y2), 
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SO 



(4.33) 



sup Vg,,(t)(?/2,00)I^^( 



> 

0<t<T 



where the third hne follows from the definition of C„ , the fact that 2/2 < 2/i < 
yl imphes Hk{y2) > Hk{yi), and the inequality F'^}{t) > - y/nt (following 
immediately from the definition of the frontier) resulting in 

(4.34) yl<y2 + V^t on {F^J (t) < ya}- 

Also, on AnnBnDCn, 



0<t<T 

= ^sup_^[vi"fcVt)(yi, 00) + pk,eHk{yi + ^/iit)]I^^(^„) ^^^^^^^ 
< -^Hk{yi) = c. 

[The second line follows from 7/2 < Ui and (4.34) and the third one from the 
definition of C^.] Thus, (4.33) yields, for n > A^, on An n Bn H Cn, 

Pkl 



and, therefore, 



A inf F'^}{t)<y2, , 



3e 



< 



^ 2(v/gfc(^ + c) ^ £ 
Pk/Hk{y2) 4 

for y2 small enough. Thus, (4.30) holds for k, j = £, y = 2/2 and n> N. 
Taking y smaller, if necessary, we extend (4.30) to k, j = i and all n. □ 

Corollary 4.7. For every j, we have 

Q'f\-oo,F^''\t)]^0, w}''\-oo,Fj''\t)]^0. 
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Moreover, for all stations j and i for which there exists a customer class k 
such that j & S{k\i), we have 

Proof. This is just a restatement of Corollaries 4.4 and 4.5, except that 
the half-line (-oo,Fj"^(t)] is now closed on the right. Corollary 4.2 asserts 

that if f]"^ were bounded below, uniformly in i G [0, T] and n, then the 
inclusion of this endpoint would make no difference. Using Lemma 4.6, we 
can ensure that Fj is bounded below with probability arbitrarily close to 
1, and the result follows. □ 

5. Inverting the frontier equations. In this section we show the first part 
of Theorem 3.2, that is, that the function <I> defined by (3.20) is a homeo- 
morphism of the set D given by (3.12) and (3.13) onto [0,cx))'' (Proposition 
5.5). It is clear that $ is continuous. Lemma 5.1 asserts that $ maps D onto 
[0, oo)*^. The proof of this lemma contains an explicit algorithm for inverting 
Lemmas 5.2 and 5.3 show that $ is one-to-one on D. Finally, Lemma 
5.4, examining the limiting behavior of $(y) in D as y ^ oo, is used to show 
that the mapping $ is open. 

Lemma 5.1. $(L>) = [0, oo)-'. 

Proof. Let w = {wi, . . . , wj) G [0, cxd)-^ be given. The aim is to find y = 
{yi, . . . ,yj) G D such that $(y) = w, that is, to solve the frontier equations 



(5.1) Wj= J2 Pk,j 
kec{j) 



Hk{yj)-HA mill y, 

iG<S(fe|i) 



j = 1,..., J, 



for a y G D. Note that if /c is a customer class entering the system at station 
j, then minig5(fc|j) yi = minig0 yi = oo, so 



Hk iVj)- Hk[ mill yi = H/, {yj ) . 

v«e<S{fc|j) 



We may rewrite (5.1) in the form 

Pk,jHkiyj)+ Y. PkJ 

keJCoU) keC{j)\Ko{j) 

For j G iTo and y G M, we define 

KoAy)= Y Pk,jHkiy)= Y PkJ 



Hk{yj)-Hk[ mill 

\ies{k\]) 



Hk{y)-Hk[ mill yi 

ieS{k\j) 
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Although defined on all of M, we shall be interested in Kqj restricted to a 
smaller set. In particular, 



^o,j-[ -oo max yl 



onto 



[0,oo) 



is strictly decreasing and has a strictly decreasing inverse 

Kq] ■ [0> oo) ^ ( -oo, max yl 

We choose ji G Jo so that K^j_^ (wj-^) = max^gj-Q KQj{wj) and we set yj^ 
^o^lil^ii)- Then 



(5.2) Wj>Koj{yj,) 
fce/Co(j) 



Vj G Jo, 



(5.3) = Koj,{yj,) = ^ Pfcji ^^fc(%i) - -f^fcf . mm m. 



Induction hypothesis. Suppose that for m = 1,...,M, we have chosen 
distinct indices ji,j2,..., Jm , have defined numbers yj^ > ?/j2 ^ ' " " ^ VjM ) 
and have defined functions 



Km-i,j{y)= X! ^fc-: 



; _ 1^ (Jl 1 ■ ■ ■ iJm — 1 ) / -x 



Hk{y)-Hk[ mill y, 

.*65(fc|j) 



ye: 



for j G J^^i'^"' Although defined on ah of M, each function 



i^m.-ij : — oo, max yl A min 



I _ 1 ■ ■ ■ )Jm — 1 ) / 



iecS{fc|j) 



onto 



;o,cx)) 



is strictly decreasing when restricted to the indicated set and, therefore, has 
a strictly decreasing inverse 



(5.4) KJ_^ . : [0, oo) ^ ( -oo, max ( A min y 



i&S{k\j) 



Suppose further that for 171 = 1, ... , M, we have G J^^iT "^""'^ and 
^m-i,i™K™)= max K-\j{wj), Vj^ = K-\jJwjJ; 
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hence, 



(5.5) 



(5.6) 



= E 



W 



pk 



fcG/C 



(il,---Jm-l)/ 

-1 u) 



Hk{yjJ-Hk[ . mill y, 



Km-l,jm iUjm) 

E 



(il.---.Jm-l)/ ■ 



(im) 



Hk{yjJ-Hk[ . mill y, 



Induction step. If M = J, we terminate the construction. If M < J, we 



proceed to step M + 1 as follows. Recall that the set J'j 



(jivJm) 



M 



is the set 



of all stations j not among ji, . . . ,jM with the property that at least one 
customer class visits j and the previous stations visited by this customer 
class are among the stations ji,-.-,jM- If there were no such station j, 
then all external arrivals would be to the set of stations {ji, . . . , jm} and all 
customers exiting a station from this set would either exit the system or else 
proceed to another station in this set. In this situation, stations outside set 
{ill • • • iJAi} would not be connected to these stations, a situation we have 
ruled out by assumption. Hence, JT^^''"'"'*^^ is nonempty. 
For j G jjf and y e M, we define 



KM,j{y) 



E 



Pkj 



k&K 



Although defined on all of 



-oo, 



(i) 



max 



4G<S(fc|j) 



U) 



yt A min y,: 



onto 



[0,cx)) 



is strictly decreasing when restricted to the indicated set and, therefore, has 
a strictly decreasing inverse 



onto 



-OO, 



max 



keK^^r-''''\j) 



We choose Jm+i G Jj 



M 



so that 



Kr/^ (wi..r.,)= max 



yi. A min y, 
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and set Vj^+i = Kj^jj^j+.C^jM+i)- Then 
Wj>KM,j{yjM+i) 



(5.7) 



J2 Pk. 



(ii Jm) 

M 



iJ) 



i65(fc|j) 



+ 



w 



Jm+1 



^MjM+iiyjM+i) 



(5i 



E 



k£K. 



(ii JA/) 



«e5{fc|jM+i) 



To complete the induction step it remains only to show that yj^^ > yj^j^^- 
We divide the analysis into two cases. 

Case I. Jm+i e jj^^lT'"''^ ■ 



In this case (5.5) implies that 
(5.9) 

fce'CAz-i ua/+i) 



Hk{yjM)-Hk[ . ^min y, 



For k e /cJ/'-'^")(jM+i) \ /c!i^':r'"-'^(jM+i), we have jm G 5(A;|jM+i), so 
(5.10) 

It follows that 



mm Hi = yjj^j . 
ie5(/c|jA/+i) 



E 



/Ofc 



(5.11) k&K^r-''^^iM\j^^:^-'''-'^UM+i) 



HkiyjM) 



— Hk ( min yi 

\i&S{k\jni+i) 



0. 



Summing (5.9) and (5.11), we obtain 

WjM+i>KM,jM+liyjM)^ 

and, hence. 



Case II. jm+i g ^if '-'^"^ \ ji^^i'''""'^- 



ACYCLIC NETWORKS 31 

In this case /C?/'-'^"^(jA./+i) / but IC^i^::'{^''-'\jM+i) = 0. This im- 
phes that for every k G /C^/' ''■''^^^(jM+i), we must have Jm £ S{k\jM+i)- 
Hence, for every k G IC^ij'"''-''^'\jM+i), equation (5.10) holds. Equation (5.8) 
becomes 

Wi.-r , , = Pk,jM+i [HkiyjM+i) - HkiVjM)]^- 



If Wjj^j_^^ > 0, then yjj^j^-^ < um- If WjM+i — 0) 'we have by definition 



= max I A mm yi 

max {yl^yjJ<yj^^, 

where the third equality follows from (5.10). The induction step is complete. 

When this construction terminates with M = J, we have chosen ji, ■ ■ ■ 
a permutation of 1, . . . , J, and we have defined numbers yj^ > — ' ' ' — Vjj > 
such that (5.6) holds for m = 1, . . . , J. Let j £ {1, . . . , J} be given, and choose 
m so that j = jm- For k G C{j) \ A^H-i '"''""^^i)) the set S{k\j) is not a subset 
of {ji, ■ . ■ ,jrn-i} and, hence, mmi(,s{k\j) Vi < Vjm = Vj- It follows that 



(5.12) E /'fcJ 

kGcu)Wi':T''"'-'\j) 



Hk{yj)-Hk\ mill w, 

i&S{k\j) 



Summing (5.6) with jm = j and (5.12), we obtain (5.1). 

It remains to show that y = (yi, . . . , yj) £ D. Let vr = (ji, . . . ,jj). By con- 
struction, for m=l,...,J, jm. G JmL'i''''"''^^ = Jm-i^ SO VT G 11. Moreover, 
Vn >yj2>---> yjj and, by (5.4), 

= ^m-l,j,„i'^jm) 



< max I A min y. 

k&K'^^^"-^''^"^'^\jm) ^ «e<S{fc|j„) 

< max vJ; . 



Thus, y G D'^ and, hence, y £ D. □ 
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Lemma 5.2. Let w = {wi, . . . , wj) G [0, oo)'^ be given. Let y = {yi, . . . , yj) 
be the solution to (5.1) constructed in Lemma 5.1 and let y = (yi, . . . ,yj) be 
another solution to (5.1). Then, for i = 1, . . . , J , we have 

(5.13) yi<yi- 

Proof. The proof proceeds by induction. Namely, let vr = {ji, ■ ■ ■ 
be the permutation constructed in the proof of Lemma 5.1. We assume that 
(5.13) holds for i = ji,. . . ,jM with some M < J (in particular, for M = 0, 
no assumption is needed). We want to show that (5.13) holds for i = ju+i- 
By (5.8) and the fact that y satisfies (5.1), we have 



^ E PKiM+i 



HkiVjM+i) - Hk[ . ^min m 



For k e ldll''"'^'"\jM+i), S{k\jM+i) ^ {ii, • • • , Jm}, so by the induction hy- 
pothesis, 

(5.15) min yj < min yi. 

«65(fc|iM+i) ie5(fc|jA/+i) 

This, together with (5.14) and the monotonicity of iJfc, yields 



Pk,j 



JM+I 



Hk{yjM+i) - Hk[ . ^min y, 



= ^MjM+liVjM+l)- 

Thus, by the monotonicity property of Kmj^j^i, either 

(5-16) yjM+i<yjM+i 

or 

(5.17) max ( y^ A min yA < y^^,^, < y^^,^, . 

keK^ll'-'"^i\j„j^,]\ »e'S(fc|jM+i) J 

However, (5.17) contradicts the definition of yj^^i^-^'- 



yjM+i = KM,ni+^ («^iA/+i ) e ( -oo, max ( yl A min y, ) 

so (5.16) holds. □ 



ACYCLIC NETWORKS 33 

Lemma 5.3. The mapping <^:D^ [0, oo)"^ is one-to-one. 

Proof. Let w = {wi, . . . ,wj) G [0, oo)-^ be given and let y = (yi, . . . , 
yj) € be a solution to (5.1). Let tt = (ji, . . . , jj) e 11 be such that y G D'" , 
in particular, 

(5.18) yh^yh^--->yh- 

In light of Lemma 5.2, it suffices to show that y and tt can be constructed 
as the output y = {yi, . . . , yj), vr = (ji, . . . of the algorithm described in 
Lemma 5.1. Once again, we proceed by induction. We assume that for some 
M < J and all m < M we have chosen in the above-mentioned algorithm 
jm = jm and j/jVn = Vjm (foi^ M = 0, nothing is assumed). We want to show 
that it is possible to choose in this algorithm yM+i and Jm+i as jjM+i and 
jM+i, respectively. By the induction hypothesis, IC^i}'"''^'^'\j) = ^\i{3) for 

ah j and J^f '-'^'''^ = jf^. In particular, ]m+i G jjf '■"'^"^ because vr G H. 
By assumption, 

(5.19) = p,-- Huitj:- ) — Hui min % 

If A; G C(jM+i)\^M '^^■'(iM+i), then5(A;|jM+i) ^ {ji, • • ■ , JA/} and, hence, 
by (5.18), min.g5(fc|j^^^^)yi <yiM+i- Thus, (5.19) reduces to 



ie5(fc|jM+i) 



^kiyi„j.,)-Hk[ min y, 



+ 



= Pk 

(5.20) 

The second equation follows from the fact that 

(5.21) S{k\jM+i) C {ji, . . . , jm} = {ji, . . . Jm} for k G /ci^^--'^")(iA/H_i) 

and, hence, by the induction hypothesis, yi = y^ for i G S{k\jM+i)- But y G 
Z)*, so 

y-- < max yt. 

By (5.18), (5.21) and the induction hypothesis, for k G )c2j'"''-'^'\jM+i), we 
have 

y-jM+i ^ y-jM = • • • = min{y,-,, . . . ,y,,J < min y,. 

*65(A;|jA/+i) 
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(5.22) 
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y- < max 



Vk ^ min Vi 

«65(fc|jM+l) 



By (5.20) and (5.22) 
(5.23) 



Since y is a solution to (5.1), for j G 



M ' 



we have 



(5.24) 



For G /C 



> E 



Hk{yj)-Hk( min y, 



Pk,j 



kGK^-'^'hj) 



Hk{yj)-Hk[ mill y, 



(ilvjA/) 

M 



(j), Siklj) C {ji, . . . Jm}, so yi = yi for z G 5(A:|i) by the 



induction hypothesis. Thus, by (5.24), 



Wj > 



E 



Hk{y-j)-Hk \ . mill y 



(5.25) 



> 



E 



The second inequahty follows from the fact that j G so j ^ 

{ill • • • )Jm} and, hence, by (5.18) and the induction hypothesis, yj < yj^^^^- 
Relation (5.25) implies 

(5.26) K^l^ito,) < K]^l^{KMAyn,J) < Vn,,,^ j ^ jt""'"^- 

By (5.23) and (5.26), Jm+i and yjj^j^-^ in the algorithm of Lemma 5.1 may 
be chosen as jpr+i and y~- . This is what we wanted to show. □ 



Lemma 5.4. liuiy^D^y^^ 



+C30. 



Proof. We argue by contradiction. Suppose that the lemma is false. 
Then there exists a sequence y„ = (y",...,yj) G D, with ||y„|| +00 as 
n ^ +00, and a finite constant M such that 



(5.27) 



|$(y„)||<M, n = l,2,.... 
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By taking a subsequence (also denoted by t/n), we may assume that for some 
vr = (ji, . . . G n, we have y„ G D'^ , n = 1, 2, .... Let 

mo = min{m G {1, . . . , J} : {yj^}n=i 2 is an unbounded sequence}. 

By definition the set D'^ is bounded above in each coordinate. Again ex- 
tracting a subsequence (still called yn), if necessary, we may assume 

(5.28) lim y" = -00. 

Let Wn = {wi, . . . , Wj) = $(yn)) n = 1, 2, . . . . By (5.6) with m = mo, we have 

(5.29) = E [h. (.J,. ) - fl. (_^ min 

Ul.---AnO-l) , 



[Recall that, by the proof of Lemma 5.3, the permutation constructed in the 
algorithm of Lemma 5.1 with input Wn = $(yn)i Vn S D'^ , may be chosen to 
be vr.] Observe that 

(5.30) S{k\jmo) ^ {ji, • • • ,Jmo-i} for k G /C^o'li^'"°~'^(jm.o), 
and by the definition of mo, 

(5.31) limsuply" I < +00, i = 1, . . . , mo — 1. 

n— >oo ' 

Relations (5.28)-(5.31) yield lim 

n^oo w'j — +00, because lim^j^—oo H}^{x) — 
+00 for all k. This contradicts (5.27). □ 

Proposition 5.5. The mapping [0,cx))'^ is a homeomorphism 

of D onto [0, cxd)"^. 

Proof. By Lemmas 5.1 and 5.3, it suffices to prove that $ is open. The 
main idea of the proof is to use the one-point (Alexandroff ) compactification 
of D and [0,oo)'' (see, e.g., [3], pages 92 and 93). Recall that the topology 
on the one point compactification X = XU {oo} of a locally compact Haus- 
dorff space X consists of open subsets of X and the complements, in X, of 
compact subsets of X. Let oo be a single point not belonging to i?'^. Let 
D = DU {oo} and M_|_ = [0, oo)"^ U {oo} be the one-point compactifications 
of D and [0,00)"^, respectively. Define ^■.D^R'^ by 



$(y), if y£D, 
00, if y = oo. 



Lemma 5.4 implies the continuity of $ at 00. Thus, <^ is a continuous map- 
ping of a compact space into a Hausdorff space and, therefore, by Corollary 
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2 on page 87 of [3], it is closed. In fact, by Lemmas 5.1 and 5.3, $ is a home- 
omorphism of D onto M^. To conclude, let U be an open subset of D and, 

hence, of T3. Therefore, ^{U) = '^{U) is open in M^. But ^{U) C [0, ooY , so 
^>([/) is open in [0,00)"^. 

□ 

6. Proofs of the main results. 



Proof of Theorem 3.2. By Proposition 5.5, only the second part of 
Theorem 3.2 needs to be shown. Let us observe that for j = 1, . . . , J, we have 
[using the convention (a, 6] = if a > 6] 



= ^(t),oo) +0(1) 



E KlmFT i^)^ min F^^it) 



+ W.7W(T^Wv.min F/")(t),oo 



+ o{l) 



(6.1) 



E 

feGC{i) L 



i(iS{k\j) 

■ E W^:^-^"^w(^^"HO,.mm 



+ 0(1) 



E/S.«)(^'"'(«),mi,.,F>"V) 



E 

kec{j) 



.min 



+ 0(1) 



ie5(fc|j) 



+ 0(1) 



ci>,-(Fr(t),...,Fj"^(t))+o(i) 



Indeed, the first equality in (6.1) holds by Corollary 4.7. The third one 
follows from the fact that class k customers with lead times at time nt 
greater than mining (^i^y-^ Fj^"'\nt) have not yet been in service at one of the 
stations i G S(k\j) and, thus, have not yet arrived at station j. Similarly, no 
class k customer with lead time at time nt greater than Fj"'\nt) has ever 
been in service at station j, so all such customers must be either in queue 
at station j or at an upstream station j'q G S{k\j). This explains the fourth 
equality in (6.1). The fifth one follows from the fact that, for every k G C{j) 
and jo G S{k\j), we have 



0<Wj;f-\t)(Fp{t),,mm Frit) 
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by Corollary 4.7. Finally, the sixth equation in (6.1) follows from Proposi- 
tion 4.1 and Lemma 4.6, together with the fact that, by definition, for every 
station i G V{k), F^'^\t) + -^/nt > y"^ and, hence, 



Hk {fI""^ {t)+^/^t)=Hk( min F^^""' {t) + ^/^t)=0. 

\tes{k\j) ) 

By Lemma 3.1, . . . , G D and, by Proposition 5.5, $ is a 

homeomorphism of D onto [0,oo)"^. Thus, (2.25), (3.21) and (6.1), together 
with the continuous mapping theorem, yield 

(F(")(t),...,Fj")(t)) = $-H(W^/")(t),...,W^J")(t))+o(l)) 

^ ^-\Wl(t\ . . . , Wm) = (FUt), • • • , F*j{t))- □ 



An), 



Proposition 6.1. Let j = l,...,J and T >0 be given. As n^oo, both 



sup sup 

iO<t<T 



n^(")(i)(y,oo)- Y: PkJ 
k&C{j) 



H,iyyFl''\t))-H,[ mm Fr it) 

J6-S(fc|j) 



and 



sup sup 

lO<t<T 



Qf (t)(y,oo)- J2 
kec{j) 

converge to zero in probability. 



Hk{yVF';''\t))-Hk{ mm Frit) 
ies{k\j) 



Proof. By an argument similar to that used to derive (6.1), we have, 
uniformly in 0<t<T, 



>v/")(t)(y,oo) 



(6.2) 



wj''\t){y\/Fj''\t),oo) + oil) 



]:^^^tyit)(yyF^-\t),mm^^^^^^ 
E 

fcGC(j) 



+ 0(1) 



vS,(0(yVi^^"^(t), min.^i^^"^(t) 



ies(k\j) 



ja&Sik\j) 



E Kj it)(y'^ Frit),. Frit) 



+ o(l) 



fcec(i) 



Y: Pk, Hk{yyFrit))-HJ mm Frit) 



+ 0(1) 



+ 0(1). 
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Moreover, the o(l) terms above may be chosen uniformly in y G M. This 
needs a justification only for the last equality in (6.2). For y > y^, the kth 
terms in the sums in both the fifth and the sixth line of (6.2) are zero. 
Proposition 4.1 gives a uniform bound for yo <y <yl., with yo arbitrary but 
fixed. Finally, the uniform bound can be extended to all y by Lemma 4.6. 
This proves the first part of Proposition 6.1; the proof of the second part is 
analogous. □ 

Proof of Theorem 3.3. Let us define a mapping ijj-.R-^ ^M'^ by 
tp{x) = {ipi{x), . . . ,ipj{x)), X = (xi, . . . ,xj), where for j = 1, . . . , J and B G 



Observe that if) is continuous. Indeed, for j = 1, . . . , J and x,y £ M , x = 
{xi,...,xj), y = {yi,...,yj), using the fact that J2k£C(j)Pk,j = 1, we have 



where A denotes symmetric difference. Therefore, by Theorem 3.2 and the 
continuous mapping theorem, we have 




sup \iPj{x){B)-^jj{y){B)\ 



Bei3{K) 




(l-Gfc(O)de 



< 2 max \xi — yi 



(6.3) 



(F,(") (t ),..., (t) ) ^ ^ (Fi* (t) FJ (t) ) . 



For j = l 



. . , J, t > and y G M 



V',-(F*(t),...,FJ(t))(y,oo) 





>V;(t)(2/,oo). 



This shows that 



(6.4) 



^(Fi*(t),...,Fj(t))=W*(t). 
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Proposition 6.1 yields 

(6.5) sup sup |W,(")(t)(y,oo)-V',(Ff)(t),...,Fj")(t))(y,oo)|40 

yeRO<t<T 

for every j = 1, . . . , J and T > 0. Combining (6.3), (6.4) and (6.5), we have 
. . . ,yvj"^) W*. The proof of {Q^^\. . . , Q^^j^) Q* is analogous. 

□ 



7. Simulation. In this section we use simulation methods to assess the 
predictive value of the theory developed in the previous sections and to 
provide a simple illustration of the methodology. In the previous sections 
we considered a sequence of queueing networks, indexed by n, whereas here 
we want to consider a single queueing network. We imagine that this single 
system is a member of such a sequence of networks corresponding to a large 
value of n, that is, a system with traffic intensities close to one. 

Here we show how the theoretical lead-time profile can be constructed 
when the system occupancy of the nth system is given for an EDF network. 
Suppressing the time variable t, we recall that we denote the queue length 

(n) 

of class k at station j by j and its scaled version by 

^(n) _ 1 ^(n) 

(n) 

We also denote the workload at station j by Wj and its scaled version by 

Recall that, by Assumption 2.1, the scaled workloads (1^/""^ . . . , Wj^^) con- 
verge weakly to (VF*, . . . , Wj). The lead-time measure Q^^^j (y, oo) represents 
the number of class k customers at station j with lead-time greater than y 
and 

Qg(y,oo) = -^Qg(V^y,oo). 

Class k customers arrive with lead-time distribution given by 

(7.1) P(4'(")<V^y) = Gfc(y). 
We define G^^\y) = Gk{^) so that 

(7.2) n4'^"^<y) = Gi"^(y) 
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is the cumulative distribution function of the lead times of class k customers 
in the nth system. The limits of the lead-time measure processes are in terms 
of the functions ff^: 

rca 

(7.3) Hk{y)^ (1-Gfc(r?))dr?. 

We also define the function 



(7.4) 



y 



Recall that the frontier at station j is and the scaled frontier is 



n 



According to Theorem 3.2, 

(7.5) (Ff \ . . . , Fj")) ^ (F;, . . . , FJ) ^ ^-\Wl,. ..,WT, 



where <^ = (<I>i, . . . 



[0, 00)"^ is defined by 



(7.6) 



fcec(i) 



Hk{yj)-Hk{ mill j/i 



j = l,...,J. 



We also define = (^>^"\ . . . , $5"^) : R-^ ^ [0, cx))-^ by 



(7.7) 



j = l,...,J. 



Suppose that, in addition to (2.17), we have 



.(n) _ 



for k E C{j), j = 1, . . . , J . Then, by (7.5) and the continuous mapping theo- 
rem, we have, for j = 1, . . . , J, 



-j=<fj (Fi ,...,F^ ) - 2^ p^^^. 



fcec(j) 



/ffc(F)-i/fc( min, F, 



= $,(F('^\...,Fr) + 0(^ 



<I>,(Fi*,...,F}) = W^; 



\3 ■ 
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Therefore, ' {Ff' F}"') « Wf' ,j = l,...,J. [The difference between 
these two quantities is 0(1), but it is small relative to the number of cus- 
tomers in the system.] Because has the same functional form as <I>, the 
proofs in Section 5 apply to as well as In particular, <I>''"^ is a home- 
omorphism of D^""^ = \/nD onto [0, 00)"^ (see Proposition 5.5). Therefore, 



(7.8) [Ft\. . . , Ff^) « (cl>(«))^^(l^f \ . . . , Wf^) ^ ( fS"\ , 

According to Theorem 3.3, for every y G M and j = 1, . . . 
(7.4), (7.5) and (7.8), 



, J, we have, by 



Q;")(y,oo)«V^Q*(^,oo 



y 



n 



(7.9) 



Hi:^\yVFp)-H, 



E A. 

fcec(i) 



min K 

ies{k\j) 



in). 



H, 



in) 



min F, 

ieS{k\j) 



in) 
■in) 



In particular. 



(7.10) 



- E Ai"^ 



in) 



min 

iG5(fc|i) 



(n) 



J = 1,...,J. 



Equations (7.9) and (7.10) indicate that the lead-time profiles can be ap- 
proximated by a deterministic function in terms of the parameters of the 
nth system, while the knowledge of the index n is not required. The above 
approximations can be verified by simulation. 



7.1. A two station case. We consider a simple network with two stations 
(J = 2) and four customer classes {K = 4). Flows 1 and 2 visit both sta- 
tions but in the opposite order, while flows 3 and 4 visit only one station 
(Figure 1). 

7.2. Constant deadline. For illustrative purposes, we present the special 
case in which customers in class k arrive at the system with constant deadline 
yl, that is, G^^\y) =Iiy*^^){y). We also assume that yl>y^>yl> y|. We 

simplify notation by writing Qj, Fj, D and Afc in place of Q^j^\ -^i"^' ^^"^^ 

in) 

and Xf^ , respectively. It is easy to see that in the case under consideration 
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k=2 
k=4 



Fig. 1. A two node acyclic network with four customer classes. 

the sets 11 and D, defined by (3.11)-(3.13), are equal to {(1, 2), (2, 1)} and 
/;){i.2) (J Di^A) ^ respectively, where 

7.11) = {{yi,y2):yi > y2,yi < yl,y2 < yl}, 

7.12) = {(yi,y2) :y2 > yi,yi < y*2,y2 < y^}- 

Given Qi and Q2, one can find Fi and F2 by inverting the system of ap- 
proximate equations (7.10), which in our case reads 

7.13) Qi « Xiivl - Fi) + X2[{y*2 - Fi)+ - (y*^ - ^2)+] + Xs{yT, - Fi)+, 

7.14) Q2 « Ai(Fi - F2)+ + A2(y2* - ^^2)+ + X^iyl - ^2)+- 

Depending on values of Qi and Q2, there are five different formulas giving 
Fi and seven for F2 presented in equations (7.15)-(7.26). These are 

Qi 



7.15) 
7.16) 
7.17) 
7.18) 
7.19) 
7.20) 
7.21) 
7.22) 



Fi^yl 



Ai' 



Fi 
Fi 
Fi 
Fi 



Ail/t + A21/2 -Q2-Q1 
Ai + A2 

^lyl + hyl - Qi 
Ai + A3 

Aiyt + A2y2 + ^3^/3 -Q1-Q2 

Ai + A2 + A3 

>^iy*i+hy*3-Qi , X2{X2y*2 + \iyl-Q2 



Ai + A2 + A3 



+ 



(Ai + A2 + A3)(A2 + A4)' 



F2^y2-^, 
A2 



F2 
F2 



Ai?/i - Q2 - Qi 
Ai 

^2y*2 + A4y| - Q2 
A2 + A4 
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(7.23) F2 

(7.24) F2 



Ai^t + \2y*2 -Q2-Q1 
X1 + X2 

Ai + A2 + A4 



F ~ ^^i^ivl + ^^yt - Qi) . ^2y*2 - Q2 
^ ' ' (Ai + A2)(Ai + A3) ^ A1 + A2 ' 

(7 26) F ~ + ^3^3 - Qi) , A2j/2 + A4I/I - Q2 

^' ' (Ai + A3)(Ai + A2 + A4) A1 + A2 + A4 ■ 

To describe the function mapping the queue lengths ((5i,Q2) to the point 
(Fi,F2) G D satisfying (7.13) and (7.14), we have divided the quadrant 
[0,00)^ in the (Qi, Q2)-plane into eight regions, I- VIII, which can be seen 
on Figure 2. Each of these eight regions is mapped onto the corresponding 
region I' -VIII' of D (i.e., I' is the image of I, etc.) plotted on Figure 3. The 
regions in [0,oo)^ are defined by the eight vertices A-H shown in Figure 2 
and the images of these vertices in D are denoted A'-H' in Figure 3. 

Table 1 gives the appropriate pair of formulas for the different ranges of 
values of (Qi,Q2), depending on the region in [0, 00)^ in which the point 
(Qi, Q2) is located. 

In the various simulation experiments, we simulate the two node queue- 
ing network, as shown in Figure 1. The external inter-arrival times and 
service times are all assumed to follow exponential distributions. For each 
simulation run, a particular queue length combination at each node, Qi = 
{Qii,Q2i,Q3i,Q4i), ^ = 1,2, is chosen. Note Qu = Q23 = 0. The simulation 
run is initiated with empty queues. Local time is accumulated when the 
queue length levels at node 1 and node 2 are exactly equal to Qi and Q2, 
respectively. At the instant the local time reaches a pre-specified value, 10 for 
the results presented in this paper, the lead-time profiles for all the customer 
at each node are recorded and the local time counter is reset to zero. The 
simulation continues until 50 lead-time profiles at each node are recorded. 



Region 


Formula for Fi 


Formula for F2 


I 


(7.15) 


(7.21) 


II 


(7.16) 


(7.20) 


III 


(7.15) 


(7.23) 


IV 


(7.15) 


(7.24) 


V 


(7.18) 


(7.20) 


VI 


(7.17) 


(7.25) 


VII 


(7.19) 


(7.22) 


VIII 


(7.17) 


(7.26) 
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Fig. 4. Profiles: Mean, Max, Mm and Theory, Qi = (50, 0, 0, 0), (52 = (20, 38, 0, 0), 
= (400,300,200, 100). 




m » TO iWf us i» m SW tm im tag w igg ISO 



Fig. 5. Profiles: Mean, Max, Mm and Theory, Qi = (40, 10, 10, 0), Q2 = (0, 30, 0, 0), 
= (200,200,110, 100). 
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The empirical lead-time profiles are expressed in the form of empirical lead- 
time c.d.f.'s. 

In the simulation we set Ai = A2 = A3 = A4 = 0.32, and iJ,k,j = 1 for all 
k,j, so that the total traffic intensity at each node is equal to 0.96. We con- 
sider three cases with different combination of queue lengths and end-to-end 
deadlines. The end-to-end constant deadline for the four customer classes 
are denoted hy D = {Di, D2, D^, D4). In Figures 4-6, the left-most dots in- 
dicate the pointwise minimum empirical cumulative distribution function of 
the lead-time profile for these 50 samples, the right-most dots indicate the 
pointwise maximum, and the central dots are the average. As a function of 
y, the ratio (Qj"^ — Qj^\y,oo))/Q^j^^ [where Q^"'\y, 00) is approximated by 
the right hand-side of (7.9)] is plotted as a dashed curve in these figures. We 
obtained the solid curves in Figures 4-6, by replacing Aj in the right-hand 
sides of (7.9) and (7.10) by Aj/0.96. This normalization by the total traf- 
fic intensity causes the theory to have better predictive value. Indeed, with 
this normalization the theoretical cumulative distribution functions and the 
pointwise average empirical cumulative distribution functions are in almost 
perfect agreement. 

In Figures 4-6, the choice of Qi and Q2 is made to illustrate different 
profile compositions at the two nodes. Figure 4 shows the profiles of the 
case when 

Qi = (50,0,0,0), g2 = (20,38,0,0) 

and 

L> = (400,300,200,100). 

In this case only flow 1 is present at node 1, while only flows 1 and 2 are 
present at node 2. (^1,^2) is solved by (7.15) and (7.23) (region III). 
In Figure 5 the queue length levels are set at 

Qi = (40,10,10,0), g2 = (0,30, 0,0) 

and the deadlines are 

D = (200,200,110,100). 

Here, all three flows are present at node 1, while only flow 2 is present at 
node 2. In this case (^1,^2) is solved by (7.18) and (7.20) (region V). 
The final case in Figure 5 shows the profiles when 

Qi = (50,0,0,0), Q2 = (50, 0,0,0) 

and 

£' = (500,100,100,100). 

Only flow 1 is present at both nodes 1 and 2. In this case (^1,^2) is solved 
by (7.15) and (7.21) (region I). 

In each of these cases the figures show the excellent predictive accuracy 
of the theory. 
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Fig. 6. Profiles: Mean, Max, Mm and 
= (500, 100,100, 100). 
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T/ieory, Qi = (50,0,0,0), Q2 = (50,0,0,0), 
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